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Abstract. The Galactic cosmic-ray propagation 
code GALPROP is designed to make predictions 
of many kinds of data self-consistently, including 
direct cosmic-ray measurements, gamma rays and 
synchrotron radiation. In the decade since its con- 
ception it has undergone considerable development, 
which is continuing. A public version was made 
available a few years ago, supported by a website. 
We describe the new features of the current version 
which will become part of the next public release. 
Plans for future developments are also mentioned. 

Keywords: propagation, software, GALPROP 

I. Introduction 

The GALPROP code (HI has now 15 years of de- 
velopment behind it. Originally written in fortran90, a 
completely new version in C++ was produced in 2001. 
GALPROP is designed to make predictions of many 
kinds of observations: cosmic-ray (CR) direct measure- 
ments including primary and secondary nuclei, electrons 
and positrons, gamma rays, synchrotron radiation, etc. 
The important concept is that these various data are 
all related to the same astrophysical components of the 
Galaxy and hence have to be modelled self-consistently. 
GALPROP is foremost a CR propagation code, solving 
the relevant equations numerically in 2D or 3D. The 
goal is for GALPROP-based models to be as realistic 
as possible and to make use of available astronomical 
information, with a minimum of simplifying restrictions. 
A complete description of the rationale and motivation 
is given in the review [18|. More details are in lfT3l . 151 , 

he mi. 

The code solves the transport equations numerically 
on a 2D or 3D grid, accounting for nuclear spalla- 
tion processes and energy losses. The solution is time- 
dependent; normally it is continued until a steady-state is 
reached, but the explicit solution in time is also useful for 
studying the effect of discrete stochastic sources. After 
the CR propagation has been computed for all species 
including secondaries, the CR spectrum at each point in 
the Galaxy is used to compute gamma rays using gas 
surveys and a model of the interstellar radiation field, 
ensuring a consistent approach. Synchrotron emission is 
also computed using the electron and positron spectra 
and a 3D model of the Galactic magnetic field. The 
main parameters for a given GALPROP model are the 



CR primary injection spectra, the spatial distribution 
of CR sources, the size of the propagation region, the 
spatial and momentum diffusion coefficients and their 
dependencies on particle rigidity. The interstellar gas dis- 
tribution is based on observations of the neutral atomic 
and molecular gas, ionized gas, while the interstellar 
radiation field (ISRF) is based on a new calculation. 
The Galactic magnetic field model can be chosen from 
among various models taken from the current literature, 
suitably parameterized to allow fitting to the observa- 
tions. 

The GALPROP code computes a complete network 
of primary, secondary and tertiary CR production as 
described in fl2l . fl51 . iflOl . starting from input source 
abundances. The code includes cross-section measure- 
ments and energy-dependent fitting functions fl4l . 151 . 
The nuclear reaction network is built using the Nuclear 
Data Sheets. The isotopic cross-section database is built 
using the extensive T16 Los Alamos compilation of the 
cross-sections [3 1 and modern nuclear codes CEM2k and 
LAQGSM [4|. The most important isotopic production 
cross-sections (2H, 3H, 3He, Li, Be, B, Al, CI, Sc, Ti, V, 
and Mn) are calculated using our fits to major production 
channels J5J, @. Other cross-sections are calculated 
using phenomenological approximations by Webber et 
al. l20l and/or Silberberg and Tsao [11] renormalized to 
the data where they exist. 

The propagation equation is solved numerically start- 
ing at the heaviest nucleus (e.g. 64 Ni), computing all the 
resulting secondary source functions, then proceeding to 
the nuclear species with A—l. The procedure is repeated 
down to A = 1. To account for some special /3 _ - 
decay cases (e.g., 10 Be^ 10 B) the whole loop is repeated 
twice. While for many purposes the 2D cylindrically 
symmetrical option is sufficient, the 3D case is more 
realistic and will be the main option used in future. The 
full time-dependent 3D solution is required for studying 
the effect of stochastic (in space and time) sources, as 
discussed in fl5l. 



II. Management and community use 

GALPROP is a public code but is in continuous 
development by a small team. It is developed and 
maintained using a variety of open source tools. Revision 
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control uses subversion^ Configuration and installation 
management is via the gnu autotoolfl The code has 
been successfully installed under a variety of unix-like 
systems, including various Linux distributions and ver- 
sions of the Mac OSX system, using different versions 
of the gcc and intel compilers. 

A websit^l with a forum for discussion of all related 
aspects is maintained, with currently about 90 registered 
users. 

A considerable number of papers in the literature 
each year make use of GALPROP. A principal current 
application of GALPROP is as the reference model 
for the Galactic diffuse gamma-ray emission and CR 
electrons for the LAT instrument on NASA's Fermi 
Gamma-Ray Observatory 0, Q], [19|. Other mission- 
related applications are CR for PAMELA, gamma rays 
for INTEGRAL, and synchrotron predictions for Planck. 
A GALPROP-based synchrotron study can be found in 
0. 

III. New features 

Features which have been added since the last release 
and which will appear in the next release include: 

• A new calculation of the Galactic ISRF using the 
FRaNKI^ coded. 

• HEALPixH output of gamma-ray and synchrotron 
skymaps. HEALPix [2] with its uniform sky coverage, 
equal-area pixels and powerful functions (convolution, 
harmonic analysis) is a standard for radio-astronomy 
applications and WMAP, Planck etc. Therefore it is 
an advantage to have this format as well as the more 
conventional (l,b) grid. 

• MapCube output for compatibility with Fermi-LAT 
Science Tools softwar^. 

• Gamma-ray skymaps output in Galactocentric rings to 
facilitate spatial analysis. 

• Implementation of alternative treatments of pion pro- 
duction. 

• 3D modelling of the Galactic magnetic field, regular 
and random, with a range of models from the literature, 
extensible to any new model as required. 

• Calculations of synchrotron skymaps on a frequency 
grid, using both regular and random magnetic fields, and 
the CR electron and positron spectra as a function of 
position (no approximating power-law assumption). 

• Parallel support with openmp 

• Memory usage optimization. 

• New HI and CO survey data, with more precise 
assignment to Galactocentric rings. 



• More accurate line-of-sight integration for computing 
gamma-ray skymaps. 

• Anisotropic inverse Compton scattering calculations 
efficiently implemented 

IV. New release 

The last public release was in 2004, and now a new 
release is planned to coincide with this ICRC. It will be 
made available via the GALPROP website cited above. 

V. Future developments 

While GALPROP represents the current state-of-the- 
art in its field, it has definite limitations, partly historical. 
For this reason we plan a completely new development. 
Features will include an adaptive grid. In addition we 
will adopt a more modern code structure taking full 
advantage of C++ and a modular design. This will 
mean that GALPROP functions will be usable by other 
codes, and new GALPROP-like codes can be built using 
the libraries provided at the users' wishes. Suggestions 
for this development from the CR community via the 
GALPROP forum are invited. 
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